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A Lagrangian from which derive the third post-Newtonian (3PN) equations 
^ . of motion of compact binaries (neglecting the radiation reaction damping) 

is obtained. The 3PN equations of motion were computed previously by 
Blanchet and Faye in harmonic coordinates. The Lagrangian depends on the 
harmonic-coordinate positions, velocities and accelerations of the two bodies. 
At the 3PN order, the appearance of one undetermined physical parameter A 
reflects an incompleteness of the point-mass regularization used when deriving 
the equations of motion. In addition the Lagrangian involves two unphysical 
(gauge-dependent) constants r[ and r2 parametrizing some logarithmic terms. 
The expressions of the ten Noetherian conserved quantities, associated with 
the invariance of the Lagrangian under the Poincare group, are computed. 
By performing an infinitesimal "contact" transformation of the motion, we 
prove that the 3PN harmonic-coordinate Lagrangian is physically equivalent 
to the 3PN Arnowitt-Deser-Misner Hamiltonian obtained recently by Damour, 
Jaranowski and Schafer. 
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I. MOTIVATION AND RELATION TO OTHER WORKS 



The long-standing problem of the gravitational dynamics of compact bodies has become 
very important in recent years because of the need to construct accurate templates for 
detecting the gravitational waves from inspiralling compact binaries in future experiments 
like LIGO and VIRGO Concerning the two-body problem the current state of the art 

is the 3PN approximation, corresponding to the inclusion of all the relativistic corrections up 
to the order 1/c^ (where c is the velocity of light) with respect to the Newtonian acceleration. 
Up to the 2.5PN or 1/c^ approximation the equations of motion are well known, as they 
have been derived by many different methods with complete agreement on the result |]^-[17|. 
They have already been used for constructing the 2.5PN-accurate templates of inspiralling 
compact binaries [p!^-pO| . 



To the 3PN order, the problem of equations of motion has been pursued by two groups 
working independently with different methods: on one hand, Jaranowski and Schafer pT|j2^ 
and Damour, Jaranowski and Schafer [p^-|25[] employ the Arnowitt-Deser-Misner (ADM) 
Hamiltonian formulation of general relativity; on the other hand, Blanchet and Faye [P^-pP|] 
work iteratively with the Einstein field equations in harmonic coordinates. Both groups use a 
regularization based on Hadamard's concept of "partie finie" to overcome the problem of the 
infinite self-field of point-like particles. However the details are actually different; notably the 
second group developed for this problem an extended version of the Hadamard regularization 
and a theory of generalized functions p7| , p8| . Both groups found that there remains one and 
only one physical constant, ct;static in the ADM-Hamiltonian formalism ||21|-25| and A in 
the harmonic-coordinate approach [p^-p9|, that is left un-determined by the point-mass 
regularization. Furthermore, in the harmonic-coordinate approach, the equations of motion 
(obtained in Ref. [^) depend on two additional constants r[ and r'2 parametrizing some 
logarithmic terms, but these constants are not physical in the sense that they can be removed 
by a coordinate transformation. The aim of the present paper is three-fold: (i) to present 
the Lagrangian of the 3PN dynamics of the compact binary in harmonic coordinates, (ii) to 
obtain explicitly from it the ten Noetherian conserved integrals of the motion in harmonic 
coordinates, (iii) to exhibit a contact transformation of the harmonic-coordinate motion 



to some pseudo-ADM coordinates in order to compare our results p6|-p9|] with the ones 
obtained by the other group [pl|-p5|. 

Concerning (i), we find a generalized Lagrangian (i.e. depending on the positions, ve- 
locities and accelerations of the bodies) whose variation yields the conservative part of the 
3PN equations of motion in harmonic coordinates as found in Ref. P^. Our second point 
(ii) is to use the fact that the Lagrangian incorporates the ten symmetries of the Poincare 
group (notably the boost symmetry) to compute the ten integrals corresponding to the en- 
ergy, the linear and angular momenta, and the center-of-mass position. In particular, we 
find that the energy agrees with the previous result of Ref. |^|. As all these integrals will 
probably be needed in future work we choose to display them explicitly, despite the length 
of the expressions. We also give the balance equations they satisfy when the radiation reac- 
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tion effect is turned on. Finally, the result of point (iii) is that there exists a unique contact 
transformation of the harmonic-coordinate dynamical variables that changes the generalized 
Lagrangian into an ordinary Lagrangian (depending on positions and velocities) whose asso- 
ciated 3PN Hamiltonian matches exactly the one given by Damour, Jaranowski and Schafer 
P^ . This proves the complete equivalence of the results obtained from the two (rather 
different) methods followed by the two groups, and constitutes a strong support of the va- 
lidity of both methods. This equivalence has also been shown independently by the other 
group (who presents also the formulas needed for computing the conserved quantities). 
Notice that it holds if and only if the un-determined constant A in the harmonic-coordinate 
formalism and the ambiguity constant cc'static in the ADM Hamiltonian are related to each 
other by 
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a result already obtained in Ref. ||26| on the basis of the comparison of the invariant energy 
of binaries moving on circular orbits. Likely the appearance of the unknown constant A 
is not due to a real physical ambiguity, but is associated with an incompleteness of the 
point-mass regularization. It is probably related to the fact that, starting from the 3PN 
order, many separate integrals constituting the equations of motion of extended bodies 
would depend on the internal structure of the objects (e.g. their density profile), even in 
the limiting case where the radius of the objects tends to zero. Further work is needed to 
compute the precise value of A. On the other hand, the constants r[ and ri^ occuring in the 
harmonic-coordinate Lagrangian disappear from the ADM-Hamiltonian (where there are no 
logarithms), in accordance with the fact that they are pure gauge. 

The plan of this paper is as follows. In Section II, motivated by the striking equivalence 
between the (regularization- related) unknown constants A and coistatic, we discuss our method 
of point-mass regularization and contrast it to the method advocated in pT|-p5|]. Section III is 
devoted to the theoretical investigations. First we recall the theory of Noetherian conserved 
quantities in the case of a generalized Lagrangian, and next we show how to eliminate the 
accelerations in the harmonic-coordinate Lagrangian by a contact transformation at the 3PN 
order. The reader interested only in the results at the 3PN order can go directly to Section 
IV, where we present the closed-form expressions of the Lagrangian and the conserved energy, 
momenta and center of mass in harmonic coordinates, and give the result for the contact 
transformation as well as the final expressions for the Lagrangian and Hamiltonian in pseudo- 
ADM coordinates. 



II. DISCUSSION ON THE POINT-MASS REGULARIZATION 



The equivalence between the respective formalisms of [PT| -P3[] and is interesting 

because the two groups have adopted some different approaches regarding the point-mass 
regularization (chosen in both cases to be based on the Hadamard concept of "partie finie" of 
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a singular function or a divergent integral |3yj3T|). Essentially the group [pi|-p5| introduced 
systematically some "ambiguity" parameters in the ADM Hamiltonian whenever the stan- 
dard Hadamard regularization yielded inconsistent results, while the group PE|-PU| looked 
for the most general solution allowed by some basic physical requirements and following 
from a new, mathematically consistent, Hadamard-type regularization. 

More precisely, in our approach [p^-|29||, we adopted some specific variants of the 
Hadamard regularization which were devised specifically for this problem |p7|,p8|. Let F 
be a function which is singular at two isolated points yi and y2, and is smooth everywhere 
else; yi and y2 are the positions of the particles in harmonic coordinates at some given 
instant t. The Hadamard partie finie of F at the point yi, denoted (-F)i, is defined as the 
angular average over all directions of approach to yi of the finite term (zeroth order) in the 
singular expansion of the function around this point. We found that this definition yields 
a natural extension of the notion of Dirac distribution at the location of a singular point, 
that we constructed by means of the Riesz delta-function ||32|. As a result, the "partie finie 
delta-function" at the point 1, denoted Pf5i where 5i = 5(x — yi), is the linear form defined 
on the set of singular functions of the type F, that associates to any F the real number {F)i 
(see Eq. (6.9) in |27|). Using an integral notation this means that J rf^x F.PfSi = {F)i. 
(The partie finie delta-function Pf5i constitutes a mathematically well-defined version of the 
so-called "good delta function" of Infeld [0.) In our derivation of the equations of motion 
at 3PN order, this prescription is employed systematically to compute all the "compact- 
support" integrals, whose integrand is made of the product of a singular potential with 
some mass density localized on the two particle world-lines. 

By applying the latter definition to the product FG we obtain f (i^x FG.Pf6i = {FG)i, 
which permits us to give a sense to the more complicated object F.PfSi = Pf (F5i), composed 
of the product of a delta-pseudo-function with a function which is singular on its support 
(such a product being ill-defined in the standard distribution theory). Namely, Pf(F5i) is 
the linear form which associates to any function G the real number {FG)i. It is important 
to realize that Pf{F6i) ^ (-F)iPf5i in general. This is an immediate consequence of the 
so-called "non-distributivity" of the Hadamard partie finie, namely the fact that (FG^\ ^ 
(F)i(G)i for two singular functions F and G in general. As an exemple taken from [jl7|, we 
have ([/^)i = [{U)iY + 2[{U)if[{U)2f , where U = Gmi/ri + Gm2/r2 denotes the Newtonian 
potential of two particles (with ri = |x — yi| and r2 = |x — y2|). In the post-Newtonian 
iteration one can check that the functions involved become singular enough so that the non- 
distributivity plays an actual role at the 3PN order: for instance, in the example above, 
will appear in the metric coefficient goo with a factor in front, which indeed corresponds 
to the 3PN order. However, there is no problem linked with the non-distributivity in the 
equations of motion up to the 2.5PN approximation [l7|. Therefore, from the 3PN order 
(but only from that order), it is a mathematically inconsistent regularization prescription 
to assume at once that J c/'^x F.Pf6i = {F)i and Pf{FSi) = {F)iPfSi. Faced with this 
problem, the authors PT|-P^ have advocated that the breakdown of the distributivity of 
the Hadamard regularization at the 3PN order is a source of ambiguities. [Actually, in 
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their first paper, see the Appendix A in Ref. ||2T|, these authors did performed their basic 
computation using the inconsistent rule Pf(F5i) = {F)iPi6i. Later in Ref. (see the 
Appendix A there), they argued that their result was "stable" against a possible violation 
of the latter rule.] By contrast, the authors pB|-pm have accepted the special features of the 
partie finie, such as its non-distributivity, and constructed by its mean a mathematically 
consistent regularization, able to give a precise sense to all computations at the 3PN order. 

The Hadamard partie finie {F)i of a singular function involves a spherical average that 
is defined within the spatial hypersurface t = const of a global coordinate system like the 
harmonic coordinates. Clearly, this definition is incompatible with the framework of a 
relativistic field theory, and we expect at some level a violation of the Lorentz invar iance of 
the equations of motion due to this regularization. Remarkably, such a violation occurs only 
at the 3PN order; up to the 2.5PN order the equations of motion in harmonic coordinates, 
as computed using the regularization {F)i, are Lorentz-invariant |]T^. To overcome this 



problem at the 3PN order, it has been necessary to define a "Lorentzian" regularization pS 



which consists merely of applying the Hadamard partie finie within the spatial hypersurface 
orthogonal to the (Minkowskian) four-velocity of a particle. It was shown in Ref. ||2^ 



that the Lorentzian regularization adds some new terms to the 3PN equations of motion 
[computed with the standard regularization (-F)i] which are mandatory in order to maintain 
their Lorentz invariance (see for instance Eq. (5.35) in [^). The Lorentzian partie finie 
of a singular function F, denoted [F]i, enables one to define a "Lorentzian" partie finie 
delta-function PfAi, namely a linear form whose action on any F gives the real number 
[F]i. It also permits the precise definition, given by Eq. (5.11) in [^, of a model for the 
stress-energy tensor of point-particles in (post-Newtonian expansions of) general relativity. 

Besides the compact-support integrals computed before, the equations of motion contain 
many "non-compact" integrals, whose support extends up to infinity and which are divergent 
at the location of the particles. To them we assign systematically the value given by the 
Hadamard partie-finie of a divergent integral: Pf J d?yiF, see Eq. (3.1) in Furthermore, 
to any F in this class, we associate the pseudo-function PfF which by definition is the linear 
form whose action on any other G gives the real number Pf j d^x FG. Given then two 
pseudo functions their product is chosen to be the "ordinary" one PfF.PfG = Pf(FG). 

An important feature of the Hadamard partie-finie integral is that the integral of a 
gradient is not zero in general, Pf J d^x diF ^ 0, since it is equal to the sums of the 
parties finies of the surface integrals surrounding the singularities when the surface areas 
tend to zero; see Eq. (3.4) in |27]. This means that the ordinary derivative of singular 
functions shows a fundamental difference with the case of regular sources, since in this case 
the integral of a gradient is always zero (provided that the integrand decreases sufficiently 
fast at infinity). One can check that some non- vanishing integrals of a gradient start to 
appear precisely at the 3PN order. Confronted with this problem, the authors [pl| -|25[] 
have considered that this signals the presence of ambiguities at the 3PN order, notably 
because their ADM-Hamiltonian density is defined only modulo a total divergence, that 
one certainly does not want to contribute even in the case of singular sources. On the 
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other hand, the authors [p^ -p9[] have accepted this feature and introduced a new kind of 
(spatial or temporal) distributional derivative acting on the pseudo-functions of the type 
Pf-F (for instance diPfF) in order to ensure that the integral of a gradient is always zero. 
It was found E7\ that it is impossible to define a derivative which satisfies the Leibniz 
rule for the derivation of a product, i.e. di{PfFG) ^ FdiPfG + GdiPfF in general, but 
that when one replaces the Leibniz rule by the weaker rule of "integration by parts", an 
interesting mathematical structure exists. By rule of integration by parts, we refer to the 
relation J d'^x [FdiPfG + GdiPfF] = 0, for F and G arbitrary functions (see Eq. (7.2) in 



27| where we use a more appropriate bracket notation for the spatial integral). While the 



rule of integration by parts is nothing but an integrated version of the "pointwise" Leibniz 
rule, the Leibniz rule itself is a stronger requirement, which is not satisfied in general as 
there are triplets of singular functions F, G, H for which j (i^x H[FdiPiG + GdiPiF] ^ 
J (i^x Hdi{PfFG). The motivation for requiring the rule of integration by parts is that it 
is clearly valid in the case of regular fluid systems. Notably it implies that the integral of a 
gradient of any singular function of type F is zero. However, because it violates the Leibniz 
rule, the distributional derivative cannot be completely satisfying on the physical point of 
view. 

Actually two different distributional derivatives, and therefore two different regulariza- 
tions, were introduced in Ref. [0. A "particular" derivative, defined by Eq. (7.7) in p7 



was first chosen for its simplicity. The two main properties of this derivative are: that 

(i) it reduces to the ordinary derivative, i.e. diPfF = Pf{diF), whenever F is bounded 
near the singularities (in addition of being smooth everywhere else), (ii) it obeys the rule 
of integration by parts. Though the particular derivative is especially convenient to use in 
practical computations, it does not follow from some "unicity" theorem. A more interesting 
derivative, on the mathematical point of view, is the so-called "correct" derivative (we follow 
the terminology of Ref. ||29|) which does satisfy a unicity theorem. Namely, this derivative 
is obtained in Theorem 4 of ||2^ as the unique derivative satisfying the properties (i) and 

(ii) above, and, in addition, (iii) the rule of commutation of successive derivatives (Schwarz 
lemma). As it turned out, the "correct" derivative, given by Eq. (8.12) of |2^, depends on 
one arbitrary numerical constant K. (Note that both the particular and correct derivatives 
reduce to the derivative of the standard distribution theory [^| when applied on smooth 
test functions with compact support.) 

Summarizing, it is possible to construct a consistent regularization based on the 
Hadamard partie finie, thus one can give a precise meaning to any integral encountered 
in the computation, but there are several possible prescriptions associated with different 
distributional derivatives (and the Leibniz rule is not satisfied). Our strategy has been 
to perform two computations of the equations of motion, associated respectively with the 
"particular" and "correct" derivatives. Then the following was shown ||2^ . 

(I) The 3PN equations of motion, when computed by means of the Lorentzian regularization 
and the particular derivative, are in agreement with the known equations of motion up to the 
2.5PN order, have the correct test-mass limit and most importantly are Lorentz invariant 
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(in a perturbative post-Newtonian sense). 

(II) Looking for the most general solution, allowed by the regularization, for the 3PN equa- 
tions of motion to admit a conserved energy and a Lagrangian description, we find that 
they depend on two unphysical gauge-constants r[ and r'2 (associated with the appearance 
of logarithms), and on one and only one physical constant A which cannot be determined 
within the method. The equations of motion possess all the physical properties that we 
expect, but the presence of the unknown constant A is somewhat baffling, as it probably 
reflects a physical incompleteness of the regularization. 

(III) When the correct distributional derivative is used instead of the particular one, the 
equations of motion depend on K in addition to r^, r'2 and A. In this case we flnd that they 
are no longer Lorentz invariant in general, but that there is a unique value of K for which 
the Lorentz invariance is recovered: K = For this value the equations of motion have 
also all the physical properties we expect. 

(IV) The different equations of motion as obtained by means of the "particular" and "cor- 
rect" prescriptions (with K = ^ m. the second case) are physically equivalent in the sense 
that they differ from each other by an inflnitesimal change of coordinates. This satisfying 
result indicates that the distributional derivatives introduced in Ref. constitute merely 
some technical tools devoid of physical meaning. 

In the scenario (III) one may wonder why after having used the Lorentzian regularization 
deflned in Ref. one still has to adjust the constant K to a. certain value in order to get 
flnally the Lorentz invariance. The likely reason is that the distributional derivatives we use 
(the particular and correct ones) have not been deflned in a Lorentz-invariant way, as their 
distributional terms are made of the delta-pseudo-function Pf5i instead of the "Lorentzian" 
delta-pseudo-function PfAi (see Eq. (3.36) in [^). As a result, we flnd in the scenario 
(III) that although most of the terms satisfy the requirement of Lorentz invariance, notably 
the terms proportional to the combination of masses m\m2 in the acceleration of particle 
one (these terms are shown to behave correctly thanks to the Lorentzian regularization), 
there still exists a limited class of terms, proportional to m|, that do not obey the Lorentz 
invariance unless K is adjusted to the value [In the scenario (I) where there is no 

constant to adjust the latter terms behave correctly.] 

The problem of the Lorentz invariance of the equations of motion was solved in a quite 



different way by the other group [l21|-P5[]. We recall that the harmonic-coordinate equations 
of motion are manifestly Lorentz-invariant because the harmonic gauge condition preserves 
the Poincare symmetry. By contrast, the coordinate conditions associated with the ADM 
Hamiltonian formalism do not respect the Poincare group, and therefore the authors [^l| -|25[| 
had to prove that their Hamiltonian is compatible with the existence of generators in phase- 
space such that the usual Poincare algebra is satisfled. More precisely, they constructed a 
generic "ambiguous" dynamics at the 3PN order, parametrized by some unknown ambiguity 
parameters associated notably with the non-distributivity of the Hadamard partie flnie and 
to the fact that the integral of a gradient, in an ordinary sense, is not zero. They showed that 
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there were only two ambiguity parameters they denoted cukinetic and ecstatic- (Actually, in the 
first paper |21|] they considered only the ambiguity constant cukinctic and obtained the value 
t^static = |- The static ambiguity was introduced in the second paper By imposing in 

an ad hoc manner the existence of the Poincare generators for their ambiguous Hamiltonian, 
they showed [El] that the parameter ciJkinetic is fixed uniquely to the value This result was 



in fact obtained earlier by comparing their expression of the energy of circular orbits ||23 



to the expression we got by means of the explicitly Lorentz-invariant formalism described 
in the scenario (I) above. Finally, having fixed cukmetic, there still remained in the ADM- 
Hamiltonian formalism one and only one undetermined constant ujsta.tici that we shall find to 
be equivalent, in the sense of Eq. ( |1.1D , to the constant A appearing in harmonic coordinates. 
[Note that, despite the resemblance between the value = ^ in the scenario (III) and the 
result cukinetic = the constant K can be fixed to this unique value only if the sophisticated 
Lorentzian regularization is used before. Without such a regularization, several other terms 
not parametrized by K would not behave correctly under Lorentz transformations, and 
therefore no value of K could be chosen in order to restore the Lorentz invariance. In this 
sense the constant K is more "specialized" than the constant co'kinetic-] 

Finally, choosing one or the other of the two approaches advocated in Refs. [^-25] and 
p6|^9[ for the regularization is a matter of taste. In view of the equivalence of the final 
results, it is a good state of affairs that the two approaches are different conceptually and 
technically. 



III. THEORY 

A. Noetherian conserved quantities for a generalized Lagrangian 

At the IPN order, the equations of motion of two compact objects in General Relativity, 
as derived in Refs. [§|,^, can be deduced from an ordinary Lagrangian, depending on the 
positions and velocities of the bodies, which was obtained by Fichtenholz [^. At the next 
2PN order, the equations of motion in harmonic coordinates, as obtained in |p|, |TO| , p 



can 



only be deduced from a "generalized" Lagrangian, depending not only on the positions and 
velocities but also on the accelerations of the particles [^. In particular, this confirmed a 
result of Martin and Sanz [p3 that A^-body systems cannot admit an ordinary Lagrangian 



description beyond the IPN order, provided that the gauge conditions preserve the Lorentz 
invariance (as it is the case for the harmonic gauge). However, it has been shown by Damour 
and Schafer that there exists a special class of coordinates, which includes the ones 
associated with the ADM formalism, such that the Lagrangian at the 2PN order expressed 
by means of such coordinates becomes ordinary, i.e. does not depend on accelerations 
anymore. This means that we can eliminate the accelerations in the harmonic-coordinate 



Lagrangian at the 2PN order by going to the ADM coordinates |T2[. In this paper, we 
shall find that the 3PN terms in the Lagrangian in harmonic coordinates depend also on 
accelerations, and that, like at the 2PN order, these accelerations can be eliminated by 
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a suitable coordinate transformation to some "pseudo-ADM" coordinates, following the 
general method of redefinition of position variables [p!^ , p5| -|37[ . 

Strictly speaking, the dynamics of two compact bodies does not derive from a Lagrangian 
at the 3PN approximation because of the radiation reaction damping effect at the previous 
2.5PN order. When speaking of a 3PN Lagrangian or Hamiltonian, we always refer to the 
conservative part of the dynamics, which corresponds to the "even" post-Newtonian orders 
IPN, 2PN and 3PN. As we shall see, the radiation reaction effect manifests itself in the 
non-conservation at the 2.5PN approximation of the conserved quantities associated with 
the conservative 3PN dynamics [see Eqs. ( |4.7| )]. 

Let us consider a harmonic-coordinate generalized 3PN Lagrangian 



L'^~'^L[yA{t),^rA{t),s.A{t)] , (3.1) 

depending on the instantaneous positions y\{t) = yA{t) (with A = 1,2 and i = 1,2,3), 
coordinate velocities v\{t) = VA(t) = djA/dt, as well as coordinate accelerations a\(t) = 
^A{t) = d\A/dt. Our harmonic-coordinate 3PN Lagrangian is given by ( [l.l] ) below, but we 
do not need to be so specific in the present Section, where most of the results hold in fact 
for A^-body systems (A = 1, ■ ■ ■ , A^). We assume that the dependence of the Lagrangian 
( ^.11 ) upon the accelerations is linear. As a matter of fact, it is always possible to eliminate 
from a generalized post-Newtonian Lagrangian a contribution quadratic in the accelerations 
by re-writing it in the form of a so-called "double-zero" term, which does not contribute to 
the equations of motion, plus a term linear in the acceleration [|12| (this argument can be 
extended to any term polynomial in the accelerations). 

The equations of motion of the Ath body are deduced from the Lagrangian by taking 
the functional derivative defined as 



.L_3L „. (3,) 



dy\ dt \ dv\ J dt"^ \ da 



"A 



We consider first, very generally, an infinitesimal transformation of the path of the particle 
A at some instant t, i.e. 5yAif) = y^l^) ~yAif)- The corresponding variations of its velocity 
and acceleration are 5v^(t) = d6yA/dt and SaA{t) = dSvA/dt. Such a transformation of the 
motion induces a variation of the Lagrangian, namely SL = L[y^,v^, a^] — L[yA,y^A, ^a] 
which is readily found to be expressible, at the linearized order in 6yA, in the form 



^^-^ + i:^%Uo(^yy . (3.3) 

where the functional derivative 5L/5y\ is given by ( p.2| ) [it is zero "on shell", i.e. when the 
equations of motion are satisfied], and where we have introduced the total time-derivative 
of a function Q = Q[6yA, Sva] defined by 
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A 



(3.4) 



Here, p\ and q\ denote the momenta that are conjugate to the positions y\ and velocities 
v\ of the particle A respectively, that is 

= (35a) 



We now discuss the Noetherian conservation laws for generalized Lagrangians following 
Refs. ||^|n|. We know from Ref. that the 3PN equations of motion in harmonic coordi- 
nates are manifestly invariant (in a perturbative post-Newtonian sense) under the Lorentz 
and more generally the Poincare group. Thus the dynamics associated with our 3PN gener- 
alized Lagrangian ( |4.1|) should stay the same after an infinitesimal Poincare transformation 
of the dynamical variables = {ct, y^). In particular, this means that = in the case of 
arbitrary infinitesimal constant spatial translations and rotations, 5?/^ = e* and = oo^jy\ 
with Uij = —ujji. In this case Eq. ( p.3[ ) implies the conservation on-shell (all the 5L/5y\^s 
are zero) of the Noetherian linear and angular momenta given by 

P^ = ^p^, (3.6a) 

A 

J' = ^^kY.iyAPA + v\q\) . (3.6b) 

A 

Thus, dP^/dt = and dJ^dt = on shell. On the other hand, we have 6L = rdL/dt in the 
case of an infinitesimal constant time translation St = r, hence the conservation on-shell of 
the Noetherian energy from Eq. ( ^.3| ), 



E = Y1 i^>A + ^\^a) - L . (3.7) 

A 

Thus, dE/dt = 0. We shall give the explicit expressions of these Noetherian energy and 
momenta at the 3PN order in harmonic coordinates in the next Section which is devoted to 
the results [see Eqs. (|3)-(|4j)]. 

Finally, let us consider the symmetry of the Lagrangian that is associated with the 
invariance under Lorentz special transformations or boosts. Clearly, since the dynamics must 
stay the same after an infinitesimal constant Lorentz boost, the corresponding variation of 
the Lagrangian has to take essentially the form of a total time derivative. At the linearized 
order in the boost velocity W^, the transformation of the particle trajectories is given by 
Sy\ = -WH + ^W^y\v\ + 0(W''W''). There should exist a certain functional of the 
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positions, velocities and accelerations such that the 3PN Lagrangian variation reads 5L = 
W^dZydt + OCW^W^), plus some "double-zero" terms at the 3PN order (which are zero 
on-shell when applying the Noether theorem). By applying Eq. ( ^.3| ), we readily find the 
conservation on-shell of the Noetherian integral = — P^t, where P* is the linear 
momentum (|3.6aD, and where represents the center-of-mass position: 



= + E {-^A + ^ Wa< + vWa + ) ■ (3.8) 

Thus, dK'^/dt = 0, or equivalently d'^G^/dt^ = (the center-of-mass vector is conserved 
in a frame where P* = 0). The existence of the latter boost-symmetry of the Lagrangian is 
a confirmation of the Lorentz invariance of the 3PN equations of motion obtained in Ref. 
p9| . The Noetherian center-of-mass in harmonic coordinates at the 3PN order is given 



explicitly by Eq. (|45| ) below. 

The ten Noetherian quantities ( |3.6| )- (|3l^ ) have been found from our generalized La- 
grangian as some functionals of the positions, velocities and accelerations of the particles. 
However, once they have been constructed, all the accelerations they involve can be order- 
reduced using the fact that they take on-shell some definite expressions depending on the 
positions and velocities as given by the equations of motion. Our final results presented in 
Section IV.A have all been order-reduced consistently with the 3PN approximation. 



B. Elimination of acceleration-dependent terms in a Lagrangian 

We start from the harmonic coordinate system = {ct, x) and perform an infinitesi- 
mal coordinate transformation to a new coordinate system x"^, generally not obeying the 
harmonic gauge condition, of the type 



x'" = x" + e^ix) , (3.9) 

where e^{x) is a function of the spatial coordinates x as well as a (local-in-time) functional 
of the trajectories yAit) and velocities v^(t) parametrized by the coordinate time t = x^/c. 
Namely, 



£^(x,t) =£'^[x;y^(t),VA(t)] . (3.10) 

Since the accelerations in the harmonic-coordinate Lagrangian appear only at the 2PN order, 
we suppose that the coordinate transformation starts at the same level. This means that 
= O (^) and = O (^) . In particular we can check that any term in the following which 
is at least quadratic in is in fact of order O (^) and thus can be neglected in our study 
limited to the 3PN approximation. The trajectories and velocities in the new coordinates 
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x'^ = {ct', x') are some functions y^(t') and v^(t') of the new coordinate time t' = x'^/c. The 
"contact" transformation of the particle variables induced by the coordinate transformation 
( ^^ -( pTTOl) is defined by Sy\(t) = y'\{t) — y\{t) (we use the same terminology as in Ref. 
jl^). Neglecting all the terms of the order of the square of we obtain 



5y^^{t)=e\yA,t)-%\yA.t) + o(^^^ . (3.11) 

In this paper we shall construct a contact transformation 6y\, composed of 2PN and 3PN 
terms and neglecting O (^), which is issued from some infinitesimal coordinate transforma- 
tion ( |3.9| )-( [3?T0| ); however we shall not be so much interested in the coordinate transformation 
itself, in particular this means that we shall not investigate to which coordinate conditions 
it corresponds to (non-harmonic and/or ADM-type). 

If the equations satisfied by the world-lines yA{t) in some initial coordinate system derive 
from the Lagrangian L, then the equations satisfied by the new world-lines y^(t') in a new 
coordinate system will derive from the new Lagrangian L' that is such that 



L'[y'Ait),<it),^Ait)Mt)] = L[yAit),VAit),eiAit)] (3.12) 

(see e.g. Eq. (5) of Damour and Schafer |jl2[)- Since we assumed that the contact transfor- 
mation SyA depends on the velocities, the new Lagrangian necessarily depends on positions, 
velocities, accelerations and also derivatives of accelerations: hA{t) = dsLA/dt. Now the same 
computation as the one leading to Eq. ( |3.3| ) shows that, at the linearized order in 6yA, 




(3.13) 



Notice that both sides of this relation are expressed in terms of the same "dummy" variables, 
chosen to be the harmonic-coordinate ones, e.g. y^. At the end, when we obtain the new 
Lagrangian, we shall have to replace this dummy variable by the one corresponding to the 
new coordinate system, y^ = + SyA- The term with a total time-derivative is the same 
as the one found in Eq. ( p.3|) , with Q given by ( |3.4| ). As one can see, the dependence 
of the Lagrangian L' upon derivatives of accelerations comes only from this total time 
derivative. Therefore, by posing L" = L' — ^ we get a Lagrangian which is dynamically 
equivalent to the Lagrangian L' and depends like L on positions, velocities and accelerations 
only 




(3.14) 



We now show that there exists a contact transformation 6yA (actually, there exist in- 
finitely many of them), together with a redefinition of the Lagrangian by the addition of 
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a total time derivative, which ehminates all the accelerations in the Lagrangian up to the 
3PN order. In other words, the 3PN Lagrangian that will follow is ordinary, i.e. depends on 
positions and velocities only. Damour and Schafer \T2\ have already shown how to eliminate 
the accelerations at the 2PN level. We shall see how to do this at the next 3PN order, but 
in fact the method is a particular application of a general algorithm to eliminate higher- 
derivative terms in a Lagrangian 1^^. Since the contact transformation (|3.11|) is assumed 
to start at the 2PN order, i.e. Sy'^ = O (^), we must control the functional derivative ^ 



appearing in the right side of Eq. ( |3.14| ) at the relative IPN order. The standard Newtonian 
contribution is then followed by a certain IPN correction, denoted m^C^, hence 



5L 



niA 



EGniB 



AB 



+ 



AB 



o 



(3.15) 



The IPN term C\ can be straightforwardly computed from the Lagrangian ( |4.1| ). The point 
is that it does depend on accelerations, C\ = C^[ys, v^, a^], with this dependence being 



linear. The presence of accelerations in C\ is the reason why the method used in Ref. ||T2 



to deal with the problem at the 2PN order cannot be extended immediately at the 3PN 
approximation. We shall see that the method necessitates the introduction in the contact 
transformation at the 3PN order of some "counter-term" X\ described below. Now, in view 



of the term — m^a^ present in Eq. ( |3.15 ), it is clear that we will be able to remove all 
the accelerations at the 2PN order if we choose for the contact transformation the term 
T^Qa ('^6 recall that q\ is the conjugate momentum of the acceleration, = Indeed 
the only possible accelerations at the 2PN order in the Lagrangian L" would be contained 
in the combination L — J^a'^aIa^ which clearly does not depend on accelerations because 



of the linearity of the original Lagrangian L upon a\. Furthermore, as discussed in [IT2 



once we have eliminated the accelerations at the 2PN order, we are free to add to the 
contact transformation any term of the type — ^r-, where F is an arbitrary functional of 
the positions and velocities only, starting at the 2PN order. This follows immediately from 
the identity ^ = ('^A^' + ^A^f")' which shows that the further accelerations produced 
by this term are contained into the total time-derivative of F, and so can be removed from 
the original Lagrangian without changing the dynamics. However, these procedures are no 
longer valid at the 3PN order because of the accelerations in the IPN term C\ of ( ^.15 ), 
which will couple to the terms ^ q\ + as suggested before and produce some new 
accelerations. The solution of the problem is to add to the contact transformation some 
correction term that we shall find to be adjustable in a unique way so that it works. 
As a result, we look for a contact transformation of the type 



Wa = — 

niA 



9a + 



dF 



(3.16) 



where q\ is defined by Eq. ( |3.5b| ); F is a general functional of the positions and velocities, 
F = F[yA, v^], and X\ denotes some "counter" term depending on positions and velocities 
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only, X\ = X\\yB,^B\- We recall that q\ is composed of 2PN and 3PN terms, which are 
easily computed from the Lagrangian ( [4.1| ). The function F must start at the 2PN order; in 
addition we assume that it contains all possible generic terms at 3PN. Finally as explained 
above the counter term X\ is purely of order 3PN. We now replace both Eqs. ( ^.15| ) and 
into L" given by ( |3.14| ) and investigate the occurence of accelerations. Among the 
terms we recognize the combination L — c^a^'a which is free of any accelerations at the 
3PN order. We also transfer several acceleration terms into the total time-derivative of F 
as before. At last we find that the only remaining accelerations in L" are contained into the 
particular combination of terms: 



X] ( ^2 





' i dF' 







+ 01^ 



As all the terms in that combination are linear in the accelerations, we see that for any 
given function F there is a unique choice of the term X\ (for each particle) such that all 
the remaining accelerations are cancelled out, namely 



B 



dF 



(3.17) 



With the latter choice, the contact transformation ( p.l6| ), defined for any F, yields a La- 
grangian L" whose only accelerations come from (minus) the total time-derivative of F. 
Therefore, the 3PN Lagrangian L'" = L" + ^ is at once physically equivalent to L", V and 
L, and free of accelerations. Our result reads then 



L"'[yA,^A] 



, i dF 
5y\ + — + 



5y_ 



dt 



(3.18) 



Remind the large freedom we still have on the definition of L'", since we constructed it for 
any functional F of the positions and velocities at the 2PN and 3PN orders. 

In this paper we shall be able to determine uniquely the function F by the requirement 
that the Lagrangian L'" be exactly the ADM Lagrangian associated with the ADM (or 
ADM-type) Hamiltonian published by Damour, Jaranowski and Schafer [^. We shall 
not give the details of the computation since it consists merely of parametrizing the most 
general function F, constructed with the dynamical variables of the problem and having a 
compatible dimension, by means of some arbitrary constant parameters, and to show that all 
these constants are uniquely fixed by the condition of matching to the ADM Hamiltonian. 
We find indeed, in complete agreement with Ref. [^], that there is a unique set of constants 
for which this works. In particular the equivalence is possible if and only if the undetermined 



constant A appearing in the harmonic-coordinate formalism p9[ is related to the constant 
i^static of Jaranowski and Schafer by Eq. ( |1 . 1| ) . Note that the latter matching shows 
also that the logarithms In (^) and In (■^) present in the harmonic-coordinate Lagrangian 
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( [4.1|) , where r[ and r!^ denote some regularization constants, are eliminated by this contact 
transformation, in agreement with the fact proved in Ref. that the logarithms, and the 
constants r[ and r2 therein, can be gauged away. See Eq. ( |4.9| ) below for the complete 
expression of the function F. 

At last, with F now fully specified by the equivalence with we obtain the ordinary 
ADM-type Lagrangian 

given explicitly at the 3PN order by Eq. ( [4.11|) below, in which, as mentionned above, we 
shall replace the "dummy" variables used in the computation, and v\, by the real dy- 
namical variables in pseudo-ADM coordinates, Y\ and V^. The ADM momentum conjugate 
to the velocity is 

,,20, 



and the corresponding Hamiltonian follows from the ordinary Legendre transformation 



^ADM ^ p^^^^^ _ ^ADM _ ^3 21) 

A 

See Eq. ( ^.12| ) for the complete 3PN expression of this Hamiltonian (as a function of and 
Pa)- [We have checked that the second equality in ( p.20|) is true at 3PN order.] Notice that, 
strictly speaking, H^^^ is not the ADM one, as it differs from it by a shift in phase-space 



coordinates at the 3PN order which is given in Ref. ||2^. Indeed, the ADM Hamiltonian at 
the 3PN order is not ordinary, as it depends on the positions and momenta as well as on their 
derivatives |^ . But this is not a concern for our purpose, since we are interested in proving 



the equivalence between our approach p6|-p9|] and the one of pT|-p5[1, that is in finding the 
existence of a unique transformation connecting both works, in whatever coordinate systems 
the two approaches found it convenient to be. We think that the equivalence found in this 
paper and in Ref. |]25| convincingly confirms the correctness of the result. This equivalence 



is especially important in view of the different procedures adopted by the two groups to 
treat the point-mass divergencies (see Section II for a discussion). 

IV. RESULTS 

A. Conserved quantities in harmonic coordinates at the 3PN order 

We first exhibit a generalized Lagrangian from which derive the 3PN equations of motion 
of two compact objects as they were obtained in harmonic coordinates; see Eqs. (7.16) in 
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29| . The Lagrangian corresponds only to the conservative part of the equations, which 
excludes the radiation reaction term at the 2.5PN order. To compute it we proceed by 
guess-work, and find the occurence of terms depending on accelerations at the 2PN and 
3PN orders. The Lagrangian is chosen to be linear in the accelerations, and to agree at the 
2PN approximation with the Lagrangian obtained in Ref. |TT[]. The result is 
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In our notation, ^-12 = |yi — y2|, ni2 = (yi — y2)/''"i2 5 and the scalar products are written 
e.g. (ni2f2) = ni2.V2. To the terms given exphcitly above, we have to add the terms 
corresponding to the relabehng 1^2, including those which are symmetric under the 
label exchange. Notice the presence of the constant A which is the only unknown physical 
parameter in this Lagrangian, and of the two unknown gauge constants r[ and (we follow 
exactly the notation of [^). The Lagrangian presented here is not the only admissible one, 
as we can always add to it an arbitrary total time derivative (double- zero terms would make 



the Lagrangian non-linear in the accelerations). We have checked that our Lagrangian ( |4.1|) 
differs indeed from the one given by Eqs. (5.4)-(5.10) in Ref. by a total time derivative. 

Next we present the expressions of the conserved integrals of the 3PN harmonic- 
coordinate motion as constructed in Section III. A. These expressions involve only the rela- 
tivistic IPN, 2PN and 3PN terms corresponding to the conservative part of the dynamics at 
the 3PN order. The radiation reaction damping effect is added afterwards. All the quanti- 
ties we present depend only on the positions and velocities, because all accelerations therein 
have been systematically order-reduced by means of the equations of motion. The energy E 
reads 
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22 2, fru\ . 22, 




-Y'^'^'^yf TY^'^''''^^''^^^ )Kl-2 + 0(-), (4.2) 



We find that this energy is in agreement with the expression obtained in Ref. by guess- 
work starting directly from the equations of motion. The logarithms In (^) and In (^) 
take the form of a gauge transformation of the energy (see Eq. (6.16) in [^). Accordingly 
they will never enter a physical result such as the circular-orbit energy when expressed in 
terms of the orbital frequency of the circular motion (see Ref. p6|). Such is not the case of 
the constant A which does enter the invariant energy. The total linear momentum P* at the 
3PN order is given by 
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Next, the 3PN angular momentum J* is 
J' = Eijkmiylv'l 
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The last constant of the motion is the vector K^. We rather present the vector G* = P^t + K^ 
which represents the center-of-mass position and varies hnearly with time. 

i I if Gmim2 , mivl\ 



1 1 / 7, ,7, ^\ 56-^^2 7G^miml 
+ ^< ViGmim2( - -jinuVi) - -(niaVa) + Z/i 7^ h 



3miuf Gmim2 f 1, 2 1/ \ 1/ \2 

H 3 ^ -{nuVi) - -[ni2Vi){ni2V2) + -{ni2V2) 

o ?"i2 \ o 4 o 

19 2 / \ 7 2 

1 f , f 235G^mlm2 f , , , A 235^277117712 A , , , 

/ 5 3 3 

+Gm-im2 ( — (ni2Vi)^ + - (7112^1 )^(7^l2^;2) + -{ni2Vi){ni2V2Y 

5 15 1 

+— (?T'i2'y2)^ - —{ni2Vi)vl - {ni2V2)vl + -{ni2Vi){ViV2) 
iZ o 4 

1 15 

+ t(^^12W2)K^2) - ini2Vi)V2 - — (ni2W2)^^2 

4 o 

+^if~T?" + — — ^f:;^(^i2^'i)^ + l{ni2Vif{ni2V2) + :^(ni2Wi)^(ni2W2)^ 
\ lb ri2 \1d 8 lb 

I 15 1 

+ 7(7T'12t'l)(ril2t'2)^ - -7:(ni2t'2)^ - 777(^12^1) V - 7^("'12t'l)(ril2t'2)t'l 

4 ib lb 2 

II 53 3 3 5 

— ^{rii2V2yvl + —v^ + -(7712^1)^(^1^2) + -^{ni2Vi){ni2V2){viV2) + 1(^12^2)2(^1^2) 

-5^;i(^;iV2) + ^iviV2y - ^{ni2Vifvl - |(ni2'i;i)(ni27;2)^^i + ^{ni2V2yv2 
o 4 o ib 

31 2 2 15, ,2 11 A ^2^2^2/79 17 

+ Y^^lf^2 - y (^l'y2)^^2 - Y^^2 ) + I Y2^^^2^l) - y (m2^^l)(ni27;2) 



17, ,2 175 2 40, , 20 2\ G^mim^/ 7 

+y (-12.2)^ - + y (^^^2) - -vij + ^:2^ ( - 3 

29, , , , 2, ,2 101 2 40, , 139 2 
+ — (^12^1) (^127^2) + 3(^i2^^2) + — - y(^^i^^2) + -^v^ 

19G^mlml G^mlm2 f 13721 22, /ri2 
H 5 -77^ ^ m ' 



8r?2 r?2 V 12b0 3 



22 



G^miml ( 14351 



12 



1260 




O 



(4.5) 



We checked that this expression of the harmonic-coordinate center of mass is changed under 
the contact transformation into the ADM-coordinate expression which is given by Eqs. (16)- 
(22) in Ref. Notice that the energy E is the only one among these integrals of the 3PN 
motion that depends on the unknown constant A. The other integrals P*, J* and G* do not 
depend on A and therefore are entirely determined. 

The latter Noetherian quantities are no longer conserved when we take into account 
the radiation reaction effect at the 2PN order. In order to express the resulting balance 
equations in the best way, we modify all these quantities by certain terms of order 2.5PN 
and find that the right-hand-sides of the equations take the form appropriate to a radiative 
flux at infinity. We pose 

2G{m2 



E 



pi 



E H — [ni2Vi] 



5c^r^ 



12 
2^2„ 



^1 



V2f + 



12 



: AGmim2 



{Vi - V2f 



2Gmi 



5c5 
Gnii 



vlv{v2 + 2—^v{v. 



Gmi 



ri2 



- mi) 

+ 1^2 
2Gmi 



+ 1 



ri2 



1^2 



{ni2Vi){v{y''2+y'A) 



12 



^1 



2G^ml 



12 



12 



+ 1^2 



G' = 
as well as 



■ 4Gmim2 ^ 
^ + Tl ^ 



{vi - V2) - 



2 2G{mi + m2] 



ri2 



— tP*. Then, the 3PN balance equations are given by 
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(4.7a) 
(4.7b) 
(4.7c) 
(4.7d) 



where the Newtonian trace-free quadrupole moment is Qij = mi{y\y{ — ^6'^^yj 



2. 



B. Contact transformation and the ADM Hamiltonian at the 3PN order 



Our final result for the contact transformation ( ^.16 ) is as follows. The first term in 



( ^.16| ) is composed of the conjugate momentum of the acceleration and is readily obtained 
by differentiating (|4.1|) : 
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The second term in (|3.16| ) involves the function F that constitutes the only possible freedom 
to adjust in order to match the harmonic-coordinate and ADM-Hamiltonian formalisms. 
This F was uniquely determined as 
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Notice the dependence of F on the logarithms, viz 
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which is necessary in order for the contact transformation to remove the logarithms of the 
harmonic-coordinate Lagrangian ([4.1|) . This result can be checked to be in agreement with 
the coordinate transformation given by Eqs. (7.2) in Ref. ||29|. The third term in (|3.16|) 
involves a correction term, purely of order 3PN, which is defined by (|3 . 1 7|) . For this term 
we get 
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The term X\ is obtained by relabeling 1^2. With those results we obtain the ADM 
Lagrangian (|3.19| ) which is an ordinary Lagrangian, not containing any accelerations, and 



furthermore not containing any logarithms. Though the investigations in Section III.B were 
done with the harmonic-coordinate quantities taken as "dummy" variables, we must present 
here the ADM Lagrangian in terms of the variables corresponding to the motion in ADM 
coordinates. We denote them exactly like in harmonic coordinates but with upper-case 
letters, e.g. R^2 = |Yi - Y2I, N12 = (Yi - Y2)/i?i2, {N12V2) = N12.V2. 
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The corresponding ADM (or, rather, ADM-type p4|) Hamiltonian is given by the ordinary 
Legendre transformation ( p.21| ) as 
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This result is in perfect agreement with the expression obtained by Damour, Jaranowski 
and Schafer (Note that in their pubhshed result, Eq. (12) in Ref. the following 
terms are missing: 



( 55 193 \ {N^2P2?Pl , , ^ 

^ „ nil ^2 h 1 ^ 2 . 

d^rl^ V 12 48 J mima 

This is a misprint which has been corrected in an Erratum p4|.) Finally, we recall that 
the agreement works if and only if our undetermined constant A is related to their static- 
ambiguity constant ecstatic by Eq. (|1.1|), and their kinetic-ambiguity constant takes the value 
i^kinetic = This completes the proof of the equivalence of the harmonic-coordinate and 
ADM-Hamiltonian approaches to the equations of motion of compact binaries at the 3PN 
order. 



ACKNOWLEDGMENTS 

We thank Bala R. Iyer for discussions. V.C de Andrade would like to thank FAPESP- 
BRAZIL for financial support. G. Faye acknowledges financial support of the EU- Network 
HPRN-CT-2000-00137. The algebraic computations reported in this paper were performed 
with the help of the software Mathematica/Mathtensor. 



28 



REFERENCES 

[1] C. Cutler, T.A. Apostolatos, L. Bildsten, L.S. Finn, E.E. Flanagan, D. Kcnncfick, 
D.M. Markovic, A. Ori, E. Poisson, G.J. Sussman and K.S. Thorne, Phys. Rev. Lett. 
70, 2984 (1993). 

[2] C. Cutler, L.S. Finn, E. Poisson and G.J. Sussman, Phys. Rev. D47, 1511 (1993). 

[3] T. Damour, B.R. Iyer and B.S. Sathyaprakash, Phys. Rev. D57, 885 (1998). 

[4] H.A. Lorentz and J. Droste, Versl. S. Akad. Wet. Amsterdam 26, 392 and 649 (1917). 

[5] A. Einstein, L. Infeld and B. Hoffmann, Ann. Math. 39, 65 (1938). 

[6] LG. Fichtenholz, Zh. Eksp, i Teor. Fiz. 20, 233 (1950). 

[7] L. Bel, T. Damour, N. Deruelle, J. Ibanez and J. Martin, Gen. Relativ. Gravit. 13, 963 
(1981). 

[8] T. Damour and N. Deruelle, C.R. Acad. Sci. Paris 293, serie II, 537 (1981). 
[9] T. Damour and N. Deruelle, C.R. Acad. Sci. Paris 293, serie II, 877 (1981). 
[10] T. Damour, C.R. Acad. Sci. Paris 294, serie II, 1355 (1982). 

[11] T. Damour, in Gravitational Radiation, N. Deruelle and T. Piran (eds.), North-Holand 
Company (1983). 

[12] T. Damour and G. Schafer, Gen. Relativ. Gravit. 17, 879 (1985). 
[13] S.M. Kopejkin, Astron. Zh. 62, 889 (1985). 

[14] L.P. Grishchuk and S.M. Kopejkin, in Relativity in Celestial Mechanics and Astrometry, 
J. Kovalevsky and V.A. Brumberg (eds.), Reidel, Dordrecht (1986). 

[15] G. Schafer, Ann. Phys. (N.Y.) 161, 81 (1985). 

[16] G. Schafer, Gen. Rel. Grav. 18, 255 (1986). 

[17] L. Blanchet, G. Faye and B. Ponsot, Phys. Rev. D58, 124002 (1998). 

[18] L. Blanchet, T. Damour and B.R. Iyer, Phys. Rev. D 51, 5360 (1995). 

[19] CM. Will and A.G. Wiseman, Phys. Rev. D 54, 4813 (1996). 

[20] L. Blanchet, Phys. Rev. D54, 1417 (1996). 

[21] P. Jaranowski and G. Schafer, Phys. Rev. D57, 7274 (1998). 

[22] P. Jaranowski and G. Schafer, Phys. Rev. D60, 124003 (1999). 

[23] T. Damour, P. Jaranowski and G. Schafer, Phys. Rev. D62, 044024 (2000). 



29 



[24] T. Damour, P. Jaranowski and G. Schafer, Phys. Rev. D62, 021501(R) (2000); and 
Erratum to be published. 



[25 
[26 
[27; 
[28 
[29 
[30 

[31 
[32 
[33 
[34 
[35 
[36 
[37; 



T. Damour, P. Jaranowski and G. Schafer, ^r-qc 001004U| . 

L. Blanchet and G. Faye, Phys. Lett. A271, 58 (2000). 

L. Blanchet and G. Faye, J. Math. Phys. 41, 7675 (2000). 

L. Blanchet and G. Faye, submitted to J. Math. Phys., gr-qc 0006100. 

L. Blanchet and G. Faye, to appear in Phys. Rev. D, gr-qc 0007051. 

J. Hadamard, Le probleme de Cauchy et les equations aux derivees partielles lineaires 
hyperboliques, Paris: Hermann (1932). 

L. Schwartz, Theorie des distributions, Paris: Hermann (1978). 
M. Riesz, Acta Mathematica 81, 1 (1949). 

L. Infeld and J. Plebanski, Motion and Relativity, Pergamon, London (1960). 
J. Martin and J.L. Sanz, J. Math. Phys. 20, 26 (1979). 
G. Schafer, Phys. Lett. lOOA, 128 (1984). 

B.M. Barker and R.F. O'Connell, Phys. Rev. D29, 2721 (1984). 
T. Damour and G. Schafer, J. Math. Phys. 32, 127 (1991). 



30 



